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Using a general Ginzburg-Landau effective Lagrangian, we study the topological structure and 
low-lying collective modes of dense QCD having both chiral and diquark condensates, for two and 
three massless flavors. As we found earlier, the QCD axial anomaly acts as an external field applied 
to the chiral condensate in a color superconductor and, as a new critical point emerges, leads to a 
crossover between the broken chiral symmetry and color superconducting phases. At intermediate 
densities where both chiral and diquark condensates are present, we derive a generalized Gell-Mann- 
Oakes-Renner relation between the masses of pseudoscalar bosons and the magnitude of the chiral 
and diquark-condensates. We show explicitly the continuity of the ordinary pion at low densities to 
' a generalized pion at high densities. 
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I. INTRODUCTION 



o ■ 

O: 

Quantum chromodynamics at finite temperature, T, and chemical potential, /x, has a rich phase structure: at low 
T and fi, the Nambu-Goldstone (NG) phase with nearly massless pions is realized by the dynamical breaking of 
chiral symmetry through condensation of quark— anti-quark pairs [l|, |2|. On the other hand, at low T and high fj,, a 
Fermi liquid of deconfined quarks is expected to appear as a consequence of asymptotic freedom Q. Furthermore, 
in cold quark matter, condensation of quark— quark pairs 0, Q leads to a color superconducting phase (CSC) and 
dynamical breaking of color gauge symmetry. At high T for any fi, the condensates melt away and a quark-gluon 
i— — i , plasma symmetric under SU(3)c is realized [6] . The experimental exploration of quark-gluon plasmas is being actively 
carried out in ultrarelativistic heavy ion collisions at RHIC, and will be continued in the future at the LHC. The 
transition from the NG phase to the CSC is also relevant to heavy ion collisions to be realized in the future at moderate 
energies at GSI, and is of interest in the interiors of neutron stars and possible quark stars. 
' We focus our attention here on the QCD phase structure at intermediate densities and study the important interplay 
■ between two competing phenomena: quark— anti-quark pairing, characterized by a chiral condensate $ ~ (??), and 
quark— quark pairing, characterized by a diquark condensate d ~ (qq). This interplay is interesting not only in its own 
, right, but also in relation to similar phenomena in other systems, e.g., the interplay between magnetically ordered 
phases and metallic superconductivity Q • Many works on the phase structure in the intermediate density region have 
J"""-. | been carried out using effective models such as that of Nambu and Jona-Lasinio (NJL) 8] and the random matrix 
model @. In Ref. [lOl ] we introduced a model-independent approach to construct a general Ginzburg-Landau (GL) 
effective Lagrangian in terms of the two order parameter fields, $ and d, and showed, in particular, that the $-<i 
coupling induced by the axial anomaly leads to a new critical point in the QCD phase diagram, implying a crossover 
between the NG and CSC phases. The crossover is relevant to the question of "hadron-quark continuity" [111 ]. 

The purposes of this paper are twofold: Firstly, we give a detailed account of the role of the axial anomaly in 
creating the topological structure of the QCD phase diagram discussed in Ref. [l(| • Secondly, we study the properties 
of the collective modes, the pseudoscalar nonet (it, K, r\, r)') and the H (the NG boson associated with the breaking 
of U(l) baryon symmetry), in the intermediate density regime where both chiral and diquark condensates are present. 

This paper is organized as follows. In Sec. |TT1 we construct a general GL free-energy in terms of the chiral and 
diquark fields that is invariant under global chiral and local color symmetry. We then use the GL free energy to classify 
possible phases for two characteristic cases: three-flavor quark matter with equal numbers of massless up, down and 
strange quarks, and two-flavor quark matter with equal numbers of massless up and down quarks; and then conjecture 
the phase diagram for realistic quark masses (light up and down quarks and a medium-heavy strange quark) from the 
results of the two cases. In Sec. IIII1 after showing the essential features of the mixing between independent massless 
modes by a toy model, we build an effective Lagrangian for the "pions" (tt, K, i]), and the H boson, and derive a 
generalized Gell-Mann-Oakes-Renner (GOR) relation in dense QCD between the masses of the pseudo-scalar bosons 
and the magnitude of the chiral and diquark-condensates. We discuss the continuity of the NG phase and CSC phase 
in the excited states. Section ITVl is devoted to summary and conclusion. In appendices, we give the detailed derivation 
of the phase diagrams, magnitudes and signs of the GL free-energy parameters, and the mass spectra of rj mesons. 
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II. PHASE STRUCTURE 

In this section, we construct, on the basis of QCD symmetry, the most general GL free-energy for the chiral field, the 
diquark fields and their mutual couplings in three spatial dimensions. We then classify the possible phase structures 
of massless QCD with two and three flavors, and discuss the QCD phase structure in the real world. 

A. Ginzburg-Landau free energy 

We write the general GL free energy, measured with respect to that of the normal phase, $ = cLl = d R = 0, as 

d L , d R ) = r» x ($) + n d (d L , d R ) + d L ,d R ), (i) 

where $, the chiral condensate, and d L R , the diquark condensates, are the basic order parameter fields; £l x d denotes 
the coupling between these order parameters. The form of O (<I>, g?l, d R ) is severely constrained by the QCD symmetry 1 

Q = SU{3) L x SU(3) R x U(1) B x U(1) A x SU{3) C . (2) 

The axial anomaly breaks the U{1)a symmetry down to Z(2Nf)A = Z(6)a where Nf is the number of massless 
flavors. Under Q, the left and right handed quarks transform as 

q L - e~ iaA e" MB V L V c q L , Qr - e iaA e~ MB V R V c q R , (3) 

where Vl,r,c are SU(3)l jR _c rotations with continuous parameters and the phases cca,b are associated with U(1)a.b 
rotations. 

We introduce the chiral fields $y (with i and j flavor indices) of mass dimension one and singlets under SU(3)c 
and U(1)b, as 

^H = G x (^ ia q i La ). (4) 

The precise proportionality constant, G x , enters only in numerical evaluation of the GL free energy parameters; see 
Appendix [B] The chiral fields transform under Q as 

$ -> e~ 2taA V L $ Vl. (5) 

In addition, 

$$t -> V L (<f><f> f )vl, -> V^($ t $)V^, det $ -> e _6iQA det (6) 

Note the invariances of $ under Z(2)a C U(1)a, and det $ under Z(6)a C J7(1)a- 

The most general form of the free energy of the chiral field, up to 0(^ A ), that is consistent with these transformations 
reads [Hi ]. 

n x = ^Tr $t$ + ^ (Tr $t$) 2 + ^ Tr ($t$) 2 _ £o ( det $ + det $t) , (7) 

where "Tr" and "det" are taken over the flavor indices, i and j. The first three terms on the right side of Eq. ((7|) are 
invariant under the full symmetry group G, while the fourth term, caused by the axial anomaly, breaks the ^(l)^ 
symmetry down to Z(6)a- The free energy tt x is bounded below for large values of the condensates for b\ + b 2 /3 > 
and &2 > 0. If these conditions are not met, we need to introduce terms C($ 6 ) to stabilize the free energy, a situation 
we will encounter later. We assume cq to be positive so that the condensate at low temperature is positive, while 
simultaneously the r]' mass obeys to 2 , > 0. We assume that ao changes sign at a certain temperature to drive the 
chiral phase transition. 



1 More precisely, to avoid double counting the Z(3) center in 5(7(3) and the discrete subgroup of 1/(1), the chiral SU(3)l,r's hi Eq. J2} 
are the quotient groups, SU(3)l^r/Z(3). 
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We define the Lorentz-scalar diquark order parameters, 

G d(q 3 Lb C Qlc) = eabceijk[<fi L ]ai, Gd(q 3 Rb C q k Rc ) = e abc e ljk [d 1 R ] ai , (8) 

where i, j, k are the flavor and a, b, c the color indices; C = ij 2 j° is the charge conjugation matrix with the properties, 
C0 = —C 2 = 1 and 7°C^ 7 = C. The 3x3 matrix [d L R ]i a belongs to the fundamental representation of SU(3)c X 
SU('3)l,r- The positive proportionality constant Gd, is introduced to make the mass dimension of d^^ R one; it enters 
only in numerical evaluation of the GL free energy parameters (see Appendix [Cj . The transformation properties of 
the d's under Q are 

d L -> e 2taA e 2laB V L d L Vj, d R -> e~ 2laA e 2ia * V R d R Vj. (9) 

Note that d^ R ) is invariant under SU(3) R ^, and Z{2) b _a(b+A)- Also the color-singlet combinations of the diquark 
fields transform as 

d L d[ V L (d L d[) Vl d R d R -» V R (d R d R ) V w (10) 

d L d R e 4M ^ (d L 4) d fl 4 -> e- 41 ^ V fl (d fl 4) V L f , (11) 

det d L -f e 6iaA e 6iaB det d L , det d fl : -► e - 6iQA e 6iQB det d«. (12) 

Note that in addition to invariances under SU(3)c, U(1)a and [7(1) B , the products d^d^ and d R d R are invariant 
under SU(3) R and 5[/(3)l respectively. Furthermore, in addition to invariances under SU(3)c, the products d^d^j 
and d R d L are invariant under U(1)b and Z(4)a C £7(1)^, while det d^ and det d# are invariant under SU(3)c, 
Z{6) B C U(1) b , and Z(6) A C U(1) A . 

The most general form of the GL free energy of the d-fields invariant under Q up to 0(d 4 ) is [HI, 0, EB|: 



0, = a Q Tr[d L d[ + d R dJj] + P 1 ([Tr(d L d[)] 2 + [Tr(d fl d^] 2 ) 

+ ft fTr[(d L d[) 2 ] + Tr[(d R d R ) 2 }) + ft Tr[(^4)(^4)] + & Tr(d L d[)Tr(d i? d^). (13) 



The transition from the normal state to color superconductivity is driven by «o changing sign. Unlike for <I>, terms 
such as det d B ,R are not allowed, since dL jR carries baryon number. 

Finally, the interaction free-energy of the chiral and diquark fields to fourth order reads [l(| 

n xd = 7l Tr[(d R d{)<l> + (d £ 4)*t] + At Tr[(d L d[)$$t + (d fl 4)*t$] 

+A 2 Tr[d L d[ + d R d R ] • Tr[$t$] + A 3 (det* • Tr^d^)*- 1 ] + fr.c) . (14) 



The triple boson coupling ~ 71, which breaks the U(1)a symmetry down to Z(6)a, originates from the axial 
anomaly. The remaining terms are fully invariant under Q. The A3 term is equivalent to the polynomial struc- 
ture, ei3k£i>j>k>®ii>$jj<(d L d R ) kkl [II]. 

Equations (O, (fl~3|) and (fT4|) constitute the most general form of the GL free energy under the conditions that 
the phase transition is not strongly first order (i.e., the values of &,d BtR are sufficiently smaller than those at zero 
temperature) and that the condensed phases are spatially homogeneous. 

B. Three massless flavors 

We first study three massless flavors. To proceed analytically, we restrict ourselves to maximally symmetric con- 
densates, namely, a flavor symmetric chiral condensate in which quark— anti-quark pairing takes place only within the 
same flavor and a color-flavor-locked (CFL) diquark condensate in which the quark-quark pairing takes place only in 
the different flavors: 

$ = diag(a, a, a), (15) 
di = -d R = diag(d, d, d), (16) 



where a and d are assumed to be real. We have chosen the relative sign between d^ and d R in Eq. (fTr?)) corresponding 
to the ground state having positive parity, as is indeed favored by the axial anomaly together with finite quark masses 
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fl7t | . The CFL condensate in Eq. (|16p is indeed realized at asymptotically high density, as shown in [l3| via the GL 
approach, and is the simplest ansatz at intermediate density. Thus the reduced GL free energy becomes 




Some comments are in order on the coefficients of the terms in ^3^(17, d): 

(i) Changes in the magnitude of a and a drive the phase transitions. In the following, we use these variables as the 
key parameters in classifying the phase structure. 

(ii) As suggested in effective theories and in lattice QCD simulations [18 ]. b may change sign as a function of T and 
fi (see Eq. (|C13|) in Appendix [C|l . To incorporate such a possibility, we have introduced a 0(a 6 ) term with a 
coefficient / > to stabilize the system. As we will see later, the qualitative phase structure for three massless 
flavors is independent of the sign of b. 

(iii) We assume that (3 is always positive, as expected from effective theories and weak coupling QCD (see Eq. (|Clip 
in Appcidix[C|- This sign implies that the transition from the normal phase to CSC is second order at the 
tree-level if there is no coupling between d and <&. Fluctuations of the diquark or gluon fields can make the 
transition first order (l9j . an effect we will not treat in this paper. 

(iv) Because both the d 2 a and tr 3 terms in Eq. (|17p originate from the axial anomaly, their coefficients, 7 and 
c, are related microscopically. Indeed, it can be shown from the instanton-induced six-fermion interaction, 
~ detij{q° R q l L ), that 7 has the same sign and the same order of magnitude as c (Appendix IB"]) . Since c is 
positive, as noted below Eq. (J7J), 7 is also positive. By lowering the free energy, the 7-term favors coexistence: 
a ^ and d ^ 0, 

(v) Microscopic calculations in weak coupling QCD and in the NJL model (Appendix [C) show that the coefficient A 
is always positive. Physically a non-vanishing a plays the role of an effective mass for the quark field, reducing 
the density of states at the Fermi surface, and the pairing energy [IH, an effect represented by \d 2 a 2 > 0. 
Furthermore, as shown in Appendix [Cl A//3 ~ ln(A/T t ;)/(/i/T ( j) 2 , which is rather small for reasonable values of 
fi, Td (the critical temperature of the color superconductivity without the a-d coupling), and A (~ /i for weak 
coupling QCD, and in the NJL model, of order the spatial momentum cutoff). Therefore, we focus primarily 
on the case 7 > with A = for the three-flavor case. The effect of the small positive A for three-flavor case is 
discussed in Appendix IA1 

In principle, the system can have four possible phases: 



The symmetry breaking pattern of the individual phases are the following: in the CSC phase, Q — ► SU(3)c+l+r x 
Z{2); in the NG phase, Q -► SU(3) C x SU(3) L+R x U(1) B , and in the COE phase, Q -> SU{3) c +l+r X Z{2). Note 
that the CFL phase and COE phase with the diquark condensate of Eq. (fTB]) each break the original U(1)b x Z(6)a 
symmetry with axial anomaly down to Z(2) (the simultaneous reflection, qL,R —> ~<1l,r), and thus these two phases 
are not distinguished by symmetry. 



The phase diagram in the a-a plane can be determined uniquely by comparing the global minima of the free energies, 
fi( NOR ), fi(° sc ), 0( NG \ and r2( COE ). In the following, we outline the basic aspects of the phase dia gram obtained 
from comparison of the free energies, and give the calculations of the phase boundaries and critical points in detail in 
Appendix 

In the absence of the a-d coupling (7 = A = 0), the four phases in (fT5]l are separated by a = (a line of second 
order transitions) and by a = a x (a line of first order phase transitions), shown in Fig[TJ where the first and second 
order phase boundaries are drawn as double and single lines, respectively. 



Normal (NOR) phase 
CSC phase 
NG phase 
Coexistence (COE) phase 



cr = 0,d = 0, 
a = 0,d^ 0, 
cr ^ 0,d = 0, 
cr ^ 0,d^ 0. 



(18) 



1. Outline of the phase diagram 
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FIG. 1: Phase structure in three and two-flavor systems without coupling between the chiral and diquark condensates (7 = 
A = 0). Phase boundaries with a second order transition are denoted by a single line, and with a first order transition by a 
double line. 

An attractive coupling, 7 > (with A = 0) leads to several important modifications, summarized in (i)-(iv) below 
and shown in FigJ5] (the left panel for small 7 and the right panel for large 7). 2 

(i) The area of the COE phase grows when both a and d 2 are non-vanishing and positive, since the —jd 2 a term 
lowers the free energy. 

(ii) The first order line between the COE and CSC phases, originally the double vertical line in Fig. [IJ terminates 
at a critical point, A, as shown in Fig. [21 This behavior is anticipated, since —jd 2 acts as an external field for a 
and washes out the first order phase transition for sufficiently large 7 or d. In Fig. [2] we denote the coexistence 
region contacting the NG phase as the NG-like COE phase, and the region which was originally the CSC phase 
as the CSC-like COE phase. 

(iii) The second order phase boundary originally located at a — splits in two, a line going to the right from the 
critical end point B and a line going to the left from the point C. Since a changes discontinuously across the first 
order phase boundary CB, the d 2 a term, which acts as a mass term for d 2 , leads to different critical temperatures 
for diquark condensation in the NG-like COE and CSC-like COE phases. 

(iv) For large enough 7, a tricritical point D appears on the boundary between the NG and COE phases. Then the 
point C, otherwise a critical end point, becomes a triple point. 



C. Massless two- flavor case 

We turn now to the massless two- flavor system (with infinite strange quark mass). In this case, all chiral and 
diquark condensates with a strange quark are suppressed. We write 

$ = diag(a,a,0), (19) 
d L = -d K = diag(0,0,d). (20) 



2 In describing the phase diagram we denote a point where a first order line turns into a crossover as a critical point; a point where a 
first order line turns into a second order line as a tricritical point; a point where a second order line terminates on a first order line as a 
critical end point; the intersection of three first order lines as a triple point; and a point where two or four second order lines meet as a 
bicritical or tetracritical point respectively [2dl |. 
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FIG. 2: Phase structure in the three-flavor system with 7 > and A = 0. Phase boundaries with a second order transition 
are denoted by a single line, and with a first order transition by a double line. On the left, C is a critical end point (0 < 7 < 
(c/3)-y^ /3/fc), and on the right, a triple point (7 > (c/3) s/ j3/b). 



The latter is the two flavor color superconductivity (2SC) state. Due to this color-flavor structure, the cubic terms in 
a and d are identically zero, and the model reduces to: 

fi 2F (a, d) = {^a 2 + V + L^j + (^d 2 + ^dA + Xd 2 a 2 . (21) 

We locate the phase boundaries and the order of the phase transitions by comparing the free energies: 
§ d 2 + fd 4 , ft( NG V,0) = fa 2 + |fj 4 + |ct 6 and fi( COE ) (a, d), with respect to ft( NOR '(0,0). 



1. Outline of the phase diagram 

We first consider b > in Eq. (|2"Tj) , with the a 6 term safely neglected, to find the qualitative phase structure. The 
system is equivalent to an anisotropic antiferromagnet (e.g., GdAlOa) [20]. For A = 0, the boundaries of the four 
phases in (jT8j) are characterized by second order lines at a — and a = with a tetracritical point at a = a = 0. In 
the presence of the repulsive d 2 a 2 term (A > 0), the area of the coexistence phase decreases, as shown in the left panel 
of Fig. [H For A > ■ky/bp, the coexistence region disappears altogether, and a first order interface between CSC and 
NG appears at a = a^f (3/b. In addition, a = a = becomes a bicritical point, as shown in the right panel of Fig. [3l 

Next we consider b < 0; here the a 6 term plays an essential role. For A = 0, the four phases in Eq. (|18|) are separated 
by a second order line at a = and a first order line at a = 36 2 /(16/) as shown in Fig. [T] With a repulsive d 2 a 2 
term (A > 0) the coexistence phase shrinks and gradually fades away as A — * 00. Moreover, a new first order line 
between NG and CSC appears and grows as A increases. This situation is shown in Fig. [4j Such a phase structure 
was previously pointed out by Vanderheyden and Jackson j21j using the random matrix model, and by Kitazawa et 
al. [22[ using the NJL model; our model-independent analysis is consistent with these results. However, the random 
matrix model does not include the physics of the suppression of the density of states at the Fermi surface leading to 
A > 0, and thus its predicted phase structure is not obviously related to that discussed here. 



D. Speculative phase structure for realistic quark masses 



The mapping of the phase diagrams obtained in the (a, a) plane to the (T, /x) plane is a dynamical question which 
we cannot address within the phenomenological GL theory. Nevertheless, on the basis of the results of Sec. [Til we can 
draw a speculative phase structure (one of many possible structures) in Fig. [5] for two light (up and down) quarks and 
a medium-heavy (strange) quark. 
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FIG. 3: Phase strucuture in the two-flavor system for 6 > and A > 0. Left: Case with a tetracritical point {^^/bp > A > 0). 
The second order line between NG and COE is characterized by a = 2a\/b, and that between CSC and COE by a — a/3/2X. 
Right: Case with a bicritical point (A > The first order line between NG and CSC is characterized by a — a\J /3/b. 




FIG. 4: Phase st rucuture of the two-flavor system with 6 < and A > 0. The second order line between NG and COE is given 
by ol = \(b - ^Jb 2 -4fa)/f. The first order line between CSC and COE is given by a = (j3/2X)[a - 3(6 - 4A 2 / ) 9) 2 /(16/)], 
while the first order line between NG and CSC is given by a 2 = (3[6bfa + (b 2 - 4/a) 3/2 - 6 3 ]/(6/ 2 ). 



In the figure, we draw two critical points: the usual one near the vertical axis (the high temperature critical point) 
originally found by Asakawa and Yazaki 18], and a new one near the horizontal axis (the low temperature critical 
point) driven by the axial anomaly. We have assumed that both critical points are in the region of positive T and /i, 
and thus the chiral transition is a crossover in the direction of both high T and high /i. With decreasing strange quark 
mass, m s , the high-T critical point approaches the vertical axis because the quark mass, which tends to weaken the 
first order transiton, becomes less effective. On the other hand, as m s increases, the low-T critical point approaches 
the horizontal axis since the system approaches the two-flavor case where the anomaly-induced critical point does 
not appear. Whether this scenario is realized or not needs to be checked in effective models incorporating the axial 
anomaly, and eventually be checked by first principles QCD simulations. Such calculations are beyond the scope of 
this paper. 

We should mention a similar two-critical-point phase structure found in the two-flavor NJL model [22]. Since the 
axial anomaly does not produce a triple boson coupling in two flavors, as discussed in Sec. Ill Cl the origin of the low-T 
critical point in Ref. [22[ is not related to ours, but can be understood as follows: Expanding the free energy of the 
two-flavor NJL model up to 0(d 2 a 4 ). we arrive at an "effective" quartic term in the GL free energy, 



-ba 4 + -bd 2 a\ 
4 4 



(22) 
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It can be shown that b changes sign from positive to negative, and b from negative to positive, as [i increases from 
zero. (See Eqs. (|C13l and lC14[) in Appendix [UP Thus for sufficiently large \x and d, the first order transition driven 
by b < disappears due to bd 2 > 0. This implies that the critical point in [22j in the two- flavor case originates from 
a large higher order term in the GL expansion, and is thus beyond the reach of our strict GL expansion. 




FIG. 5: Schematic phase structure with two light (up and down) quarks and a medium heavy (strange) quark. The arrows 
show how the critical point and the phase boundaries move as the strange-quark mass increases towards the two-flavor limit. 



III. EXCITATION SPECTRA OF MESONS 

We turn, in this section, to study the mass spectra of mesons including H boson in the intermediate density region 
for the degenerate three-flavor case. (The rj is discussed in Appendix [D]) We focus on the region near the phase 
boundaries where a and d are sufficiently small in the COE phase, and consider energy scales, p, smaller than the 
pairing gap d, so that we can neglect excitations of the quarks. 

The symmetry breaking pattern for the chiral field at low density is SU(3)l x SU(3)r — > SU(3)v- Accordingly one 
has eight Nambu-Goldstone (NG) bosons, identified with the pseudoscalar qq-state octet. These pions are massless in 
the chiral limit, with masses generated by finite quark masses. (The rf meson is much heavier than the pions because 
of the U(\)a axial anomaly.) 

On the other hand, the pattern of symmetry breaking for the diquark fields in the asymptotic high density limit is 
SU(3) C x SU(3) L x SU(3) R x U{l) B -> SU(3) c +l+r x Z(2) b , leading altogether to 17 NG bosons. Eight of these 
bosons are "eaten" by the gluons through the Anderson- Higgs mechanism; the remaining nine are the pseudoscalar 
octet 7r (the "generalized pions," or qqqq-states as discussed in [1, [IS US EL HH] ) , together with the H boson. 3 

In the intermediate density region, two sets of light pseudoscalar mesons are present, one arising from the chiral 
condensate (qq-states) and the other from the diquark condensates (qqqq-states); these two sets are mixed through 
the interaction terms [23 |. Specifically, we have two sets of pion octets (n and 7r) (as well as rf and ff) from the chiral 
and diquark fields, and in addition the H boson from the diquark field in this region. As we show below, mixing 
causes one of the two pions to become heavy while the other becomes massless. (A similar discussion of the mixing 
between qq- and qqqq-states at T = and /i = is given in [26j].) The H boson does not experience mixing, since it 
exists only in the diquark sector, but remains exactly massless independent of the quark masses. 

In the following, we examine the mass spectra corresponding to this physical picture more quantitatively using 
an effective Lagrangian approach. We begin with a toy model with independent massless fields to show explicitly 
the essential features of the role mixing plays. Next we construct an effective Lagrangian involving the pion octet 
appropriate in the intermediate density region. Using the nonlinear chiral Lagrangians studied for the NG phase 



In addition, one also has a massive fj' , since U(1)a remains broken by the axial anomaly in the high density limit. However, the violation 
of (7(1) a symmetry becomes much smaller as screening suppresses the instanton density by a high power of [17j |. Correspondingly 
we have a light rf meson at asymptotically high density. 
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[13] and for the CFL phase [U [H, [29|, we take into account ir and ff and study the effects of the mixing 
terms determined in Sec. [Ill which were previously not considered. The diagonalization of the mixing matrix for tt 
and tt leads to a generalized Gell-Mann-Oakes-Renner (GOR) relation in dense QCD, which relates the masses of 
pseudoscalar bosons to the magnitude of the chiral and diquark condensates. We also discuss the continuity between 
the NG and the CSC phases associated with the excited states. 



A. A toy model 

To show the effect the mixing of independent fields in the presence of mass terms, we first consider a simple toy 
model for two complex scalar fields, a and a: 



mix i '-'mass i 



, 

2 2 
A 

mix 



£kin = — |<V*I + T l<^<S 



C m w = — (aa^ + h.i 



771 

C mass = -(Ba + Ca + h.c), (23) 



where A, B, C, f, f and m are positive constants. 

We focus on the low mass phase excitations, described by the fields <p and (p defined by 



= exp (itp/f) , a = exp (i<p/f \ . (24) 



In the absence of the mixing term £ m i x and mass term £ mass , both ip and <p are massless. We turn on £ m ; x and £„ 
and expand the Lagrangian in terms of ip and (p. After subtracting a constant part, we obtain 



(^ + toB)// 2 -A/ ff 
-A/ff (A + mC)/p 



Mr \ — \ \^^" lD )n jj (2 K\ 



The matrix Mtu>\ is diagonalized by a transformation 



y>i \ / cos sin Q \ ( (p 
<f2 J \ — sin cos 9 J \ ip 



(26) 



where O is the mixing angle; we thus we obtain the eigenmasses mf pi and m 2 2 {mj pi < WC 2 ): 
(i) For m = 0, we have 

r<(m = 0) = 0, (27) 
= 0) = ^fi + iV (28) 



V/ 2 z 2 / 

The presence of a massless excitation corresponds to simultaneous rotation of ipi and y?2 in the same direction, while 
the massive mode corresponds to rotation in opposite directions. 



(ii) For small but non-vanishing m, we obtain 



in 



,2 



^-AB + C), (29) 



p + P 

ml ~ m 2 2 (m = 0) + + C^J , (30) 

In Sec. MI Cl we will discuss mass relations of this kind arising via mixing in the realistic situation of QCD. In particular, 
Eq. (12"9")) is a key relation corresponding to our later generalized Gell-Mann-Oakes-Renner relation, Eq. ([56]) . 
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B. General effective Lagrangian 

To derive the effective Lagrangian for the light excited states, the pions (tt, K and 77), 77' and the H boson in 
the intermediate density region, we fix, as above, the magnitude of the chiral and diquark fields and consider only 
fluctuations of their phases about their vacuum configurations. We thus parametrize the fields as 

$ = a£e^ 9 , d L =dU L e 2 « §+ *\ d R = -dU R e- 2 ^~*\ (31) 

where £, Ul and U R are SU(3) matrices, and the angles 9 and 9 are U{\)a phases, and </> the U(1)b phase. 

As mentioned at the top of Sec. IIII1 all eight gluons acquire a mass of order gf^ ~ O(gn) by "eating" the eight 
colored fluctuations of Ul.Ri where g is the QCD coupling constant and /# is the decay constant associated with 
Ul.r, defined below. On the low momentum scales we consider, p < d <C gf* , gluons are not low-lying modes. The 
remaining color-singlet fluctuations correspond to n and are parametrized by the field 

A = U L U R , (32) 

which transforms under SU(3)l X SU(3) R as 

A -» V L AVl (33) 



1. Kinetic terms 



The kinetic term invariant under Q for the chiral fields £ and 9, to second order in derivatives is, 

4" = /'ST 1 * (^SS.Et) + k- g ^d^9d y 9, (34) 



where the metric tensors g% v = diag(l, v 2 ., v 2 , v 2 ) and g^Y = diag(l,i> 2 ,,7j 2 ,,u 2 ,) arise due to the absence of Lorentz 
invariance in the medium; Here v n and ly are the speeds of the pions and 77', respectively. 
The standard pion fields tt 3 are defined by 

E = cxp (iX j n j /U) , (35) 



where the A J (J — 1, . . . , 8) are Gell-Mann matrices normalized so that TrA l A J = The first term in Eq. 

the standard leading-order chiral Lagrangian j27| , except that the speed of the mesons in the medium differs from the 
speed of light and the Lagrangian contains additional contribution from the field 9. 

Following the discussions of Refs. [1 61 ] and [24|, the most general Lagrangian of the diquark kinetic term invariant 
under Q with two derivatives is 

t%* = (flj.A&At) + ^YdJdJ + ff&d^t, (36) 

with the additional metric tensor g^Y = dia,g(l,v'j 1 ,v'j 1 ,v'j 1 ). The generalized pion fields tt 3 are defined by 

A = cxp . (37) 

The three terms in Eq. (f3"6"|) are the kinetic terms for pions, 77', and the H boson, respectively fill I24I [28j . 

At extremely high density, the decay constants fa , and fn are found by matching to their microscopic asymptotic 
values [ill ]: 

2 21 - 81n2 m 2 2 

J* lg 2 ^ 2 ' lH *' 4 27T 2 ' 1 ' 

while the velocities approach the asymptotic sound speed: t; 2 - H — > 1/3. 
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2. Mass terms 



We now assume small but finite quark masses and construct the possible terms in the GL free-energy allowed by 
the symmetry group G, where M is the 3x3 quark mass matrix in flavor space. In the intermediate density region, 
where U(1)a symmetry is violated by the axial anomaly, the lowest order possible mass terms in the GL free-energy 
are O(M), unlike at asymptotically high density where the leading terms are 0(M 2 ). 

The bare quark mass term in QCD reads 



£qcd = Q L Mq R + h.c. 



(39) 

As in the standard procedure to build up the chiral Lagrangian, we treat the matrix M as a spurion field and assume 
it to transform under Q as, 



M 



e - 2iaA V L MVl 



'r- (40) 
This transformation law together with Eqs. ([5|), (flQ|) and (TTTj) enable us to write down the possible mass terms. It 
should be stressed that $, dtd^ and M share the same transformation property under Q except for the U(1)a phase 
rotation. Therefore, the guiding principle to construct the quark mass terms is to replace the $'s or d^d^s involved 
in the GL potential Q x ,d, x d, in Eqs. 0, $F5§ and with M. 

The general mass terms of the chiral and diquark fields obtained from this procedure are: 

+ B [e abc Ei jk M ai $ bj $ ck + h.c] 
+ Ci [Tr (Af$t) Tr + h.c] 

[Tr{(M&) ($$ f )} +h.c] 



C 2 
Ti 

Ai 

A 2 

A, 



)} 



d L d\ 



+ d R d R 
d R d R 



Tr {m ( 

-{( 

Tr (c 

£abc£ijkM ai $ bj (d L d R ) + 

\ / ck 



h.c. 

)(M$t)} + 
Tr (M&) + 
h.c. 



h.c. 
h.c. 



(41) 



Note that the terms proportional to Bq and Li, which break U(1)a symmetry, originate from the axial anomaly, or 
equivalently, by instanton-induced interactions. The leading mass terms, which are linear in M and have minimal 
number of the chiral and diquark fields, become 



£ mass = Aq [ Tr ( M$ t) + h.c] + Ti Tr hd (d R dty } + h 



(42) 



3. Overall effective Lagrangian 



We turn now to the overall effective Lagrangian for 7r, rj and H. With the parametrizations of the fields defined in 
Eqs. (|3ip and (|3"2")1 . the GL potential of the chiral part Vl x in Eq. ([7]) and the interaction part Q, x d in Eq. (fT4")) reduce 
to 



£ x = cqct 3 cos (69) + const. 



£ x d = "i\d 2 o 



+ X 3 d 2 a 2 



3 2i(e-e) T r (^t s 



const. 



£ cff = 



reft 



reft 

'Hv') 



M, 



?7' 1 



(43) 
(44) 



while the diquark part f2<z in Eq. ([13]) reduces to a constant and can be neglected. From Eqs. (|34|) . (|36|) , (|42[) . (|43|) 
and (|44p . the overall effective Lagrangian takes the form up to 0(ir 2 , tt 2 , r\ 2 ,f\ 2 ): 



(45) 
(46) 

(47) 
(48) 
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FIG. 6: Mass spectra of pions in the intermediate density region: (a) without a-d mixing (7 = A = 0); (b) with a-d mixing 
(7 / 0,A / 0); and (c) in addition with finite quark masses (m q 7^ 0). Explicit expressions for the m^-'s are given in the text. 



where the quark mass matrix is assumed to be diagonal and flavor symmetric, M = diag(m ? , m q , m q ). The fields are 
redefined as 



fv' ' fw ' /h 



(49) 



The mass matrices of (tt, tt) and (77', rj') are 

/ jr(7id 2 <J + X 3 d 2 a 2 + A Q m q a) -jlj-facPa + X 3 d 2 a 2 ) \ 

Mm ~ [ ' -jjjrhid 2 * + \ 3 d 2 a 2 ) ^( 7l tpJ + A 8 (?o 2 + r 1 m,d a ) J' (50) 

■jr (3c ct 3 + 2 7 id 2 cr + 2A 3 d 2 a 2 + 2A a m q a) j-^iA^cfa - 2\ z d 2 a 2 , 

M ("'' = 12 I 1 (4 7l d 2 <7-2A 3 d 2 (T 2 ) ■±(8^ 1 d 2 a + 2\ 3 d 2 a 2 + 8T 1 m q d^ ' ' (51) 



Since we have degenerate masses of up, down, and strange quarks, the Bedaque-Schafer term [3fj | originating from 
the effective modification of the chemical potential due to the mass differences does not appear here. The off-diagonal 
components of Mi v \ originate from the chiral-diquark coupling. As the form of Mm shows, both 7r and tt are exactly 
massless in the absence of a-d coupling and quark masses. 



C. Mass spectra of low- lying collective modes 



In this subsection, we investigate the mass spectra of pions and H on the basis of the effective Lagrangian obtained 
in the last section. (For 77', see Appendix [Dl) In particular, we obtain a generalized Gell-Mann-Oakes-Renner (GOR) 
relation in the intermediate density region. We also discuss a possible connection to hadron-quark continuity for these 
excited states. 



1. H boson 



As we see from Eq. (|47|) . the H boson is exactly massless irrespective of the quark masses. This is because the H 
boson only exists in the diquark part and is not affected by the chiral-diquark coupling, so that there are no mass 
terms for H involving the quark masses. On the other hand, the mass spectra of tt and 77', unlike the H boson, are 
non-trivial on account of mixing through the chiral-diquark coupling. 



2. Pions 



The eigenstates of the mass matrix Mf w \ in Eq. ([50]) can be written as 

7Tl \ _ / COST? sint? \ / TT \ . . 

7T 2 \ - Sin 7? COST? J \ 7f I ' 1 ' 
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with mixing angle z?. We compute the pion masses m 7Vl and m n2 (we take m 2 . i < m 2 2 ) with and without quark 
masses. 

(i) In the chiral limit, m q — 0, the pion masses are 

m^m^O) = 0, (53) 
<K=0) = MV + A 3 dV)^i-A (54) 



with the mixing angle $q, which satisfies 

tan ?9n = 

/tt 



tantfo = x". (55) 



(ii) For non- vanishing m 9 , we find 

,2 



with the mixing angle 



^^(A a + T ld 2 ), (56) 
ml 2 = ml 2 (m, = 0) + ^ AqCT + j§r ld 2 ) > (57) 



tan^ = — + — — — — mq. 58) 

A /* (/£ + f£)hid 2 a + X 3 d 2 a 2 ) 



The resultant mass spectra of pions is shown in Fig. [Sj Panel (a) shows the case with neither the chiral-diquark 
coupling nor finite quark masses. In this case both tt and tt are massless and degenerate. Figure [6](b) shows the 
effects of the chiral-diquark coupling without finite quark masses; here the masses are split, and a massless mode 
still remains. Panel (c) shows the effects of both chiral-diquark coupling and finite quark masses; here because of 
the explicit breaking of chiral symmetry by quark masses, the masses of ix\ and tt-i are lifted by AmJ and Am^ 2 
respectively. The mass formula for the lightest pseudoscalar meson, Eq. (|56|) . is a generalized form in dense QCD of 
the Gell-Mann-Oakes-Renner (GOR) relation connecting the masses of pseudoscalar bosons to the chiral and diquark 
condensates. The T\ term, which originates from the effect of the axial anomaly in the pion Lagrangian, 
shows the crucial role of the axial anomaly not only for the phase structure but also for the excitation spectra in 
the intermediate density region. When the diquark condensate d decreases as the density becomes low, Eq. ([56]) 
reduces to the standard GOR relation f%m 2 = A[)m q a. On the other hand, at asymptotically high density, the chiral 
condensate a is small and the axial anomaly is highly suppressed as Ti ~ /i(Aqcd / m) 9 (1 / ff) 14 [l3- Then, the linear 
term in m q disappears on the right side of E q. (1561) and the leading term becomes 0(m 2 ), with m 2 i oc m 2 d 2 . This 



UCllll 111 l'l*q »-llOCipjJCCLl D Ull tllC llglll OIVIC Ul 1~IK±. \IU\JII dlll_l llll_. l^diaill^, 

result is consistent with observations given in Iff 31 1 at high density. 

Note that 112 is much heavier than 7Ti and it is thus expected to be radically unstable against decay into tti. The 
surviving excitation mode tt\ exhibits hadron-quark continuity with increasing baryon density. The mixing angle 
d = describes the low density limit (jti — n) and 1? = tt/2 the asymptotically high density limit (tt± = tt). In the 
intermediate density region, w± is the result of mixing between the a ~ qq state, as seen in the low density region, and 
the d 2 ~ qqqq state, as seen in the high density region. The mixing angle $ changes with the chemical potential fi. 
Therefore, across the entire span of the chemical potential, we have light pseudoscalar modes associated with chiral 
symmetry breaking, indicating hadron-quark continuity for the excited pionic states. 



IV. SUMMARY AND CONCLUSION 



In this paper, we have investigated the phase structure and excited spectra of dense QCD in the presence of the 
axial anomaly. We constructed the general GL free-energy and effective Lagrangian for the chiral, diquark fields, and 
their interactions, consistent with QCD symmetry, Q = SU(Nf)L X SU(Nf)n x U(1)b X U(1)a x SU(S)c- 

In Sec. HH we have shown that the QCD axial anomaly acts as an external field for the chiral condensate under 
the influence of the diquark condensate. The first order chiral transition is changed to crossover for a large diquark 
condensate, and a new critical point driven by the axial anomaly emerges in the QCD phase diagram. The former is a 
realization of hadron-quark continuity for the QCD ground states. Determining the precise location of this new critical 
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point is a future task for phenomenological models and lattice QCD simulations. Our schematic phase diagram would 
be made more realistic by including effects such as finite quark masses, charge neutrality, /^-equilibrium, and thermal 
gluon fluctuations (TEL . Open questions include whether the new critical point would survive in an inhomogeneous 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, and how the COE phase at low T and low fi, Fig. [5j is affected by 
quark confinement. 

In Scc- IIII) we have shown that the massless pions in the chiral limit in dense QCD result from the mixing of it ~ qq 
and 7r ~ qqqq states. Turning on a small but finite quark mass, we found that the pions acquire a mass both from the 
chiral condensate a and the diquark condensate d 2 , which is summarized in the generalized Gell-Mann— Oakes-Renner 
formula, Eq. (|56|) . Here the QCD axial anomaly plays a key role, giving the ^-contribution to the pion mass. Our 
results show explicitly that hadron-quark continuity is realized not only in the ground state but also in the low-lying 
collective excitations. The continuity of the heavier excitations such as the vector mesons and baryons remains an 
interesting problem for future exploration. 

In Sec. IIII1 we focused on the COE phase with three degenerate flavors. However, in the intermediate density 
region, in which we are interested, flavor symmetry breaking due to m s 3> m u ,d is non-negligible and causes mass 
splittings of the octet of pions in a non-trivial way. In particular, it would be interesting to study how, across the 
intermediate density region, the inverse mass ordering, m^t < ~ < 16], at high density, turns into the 
normal mass ordering, < nix < rn v < m v ' , at low density. 
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APPENDIX A: DERIVATION OF PHASE BOUNDARIES AND CRITICAL POINTS 



In this Appendix we describe more quantitatively the phase boundaries for three and two massless flavors respec- 
tively for both a positive and negative fourth order term, b > or b < 0. We only outline the derivation of the phase 
boundaries and critical points; the results can be confirmed explicitly by comparing the free energies. 



1. Three massless flavors 



To find the minimum of the GL free energy Sl^p in Eq. (TT7|) , it is useful to eliminate either of the variables a or d 
by solving the stationarity condition, 



«^3F fa i p 2 



3d = 2d U + 2 d - 7CT J=°' (M) 

which yields d — and d 2 — (2//?)( 7 cr— a/2). The phase boundaries and the order of the phase transitions are obtained 
by comparing the free energies: fi( csc )(0, d)(= ad 2 /2 + /3d 4 /4), ft( NG ) (a, 0)(= aa 2 /2 - ca 3 /3 + ba 4 /4 + /ct 6 /6), and 
J7(coe) measurec ] w ith respect to that of the normal phase, f2( NOR ) (0, 0). Using the d ^ solution of Eq. (|A1|) . we 
have 

^ (COE V, d(a)) = -g + 7 V + ^ 2 - |o» + ^ 4 + f 6 (<r > ^) , (A2) 



with 



or equivalently, 



with 



7* = "7//?, (A3) 
a* = a-2-f 2 /(3, (A4) 



o (COE) (cr(d), d) = n + + Y rf4 + higher order terms ' ( A5 ) 



fi = (48 7 4 a - 16 7 3 ca + 6 7 2 6a 2 + /a 4 ) , (A6) 

(/a 4 + 4 7 2 6a 2 - 8 7 3 ca + 16 7 4 a) , (A7) 
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a. Positive b 



For b > 0, the sixth order / term in Eq. (fl~7| does not change the qualitative structure of the phase diagram and 
can be safely neglected. In the following, we study the phase boundaries ABCC, BB', CD, and CDD' in Fig. [2] 
separately and find the locations of the characteristic points A, B, C and D. 

• The critical point A and the first order phase boundary ABCC 

We first eliminate the er 3 term in Eq. (|A2|) for the COE phase by introducing a new field r = a — c/(36), in terms 
of which, 



^ COE )(r)=« c+7c V+|r 2 + V, 



(A9) 
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with 



" c = ^ + 36 7 + 18V a ~108p' 

+ (M0) 
° c = " 36' 



where a* and 7* are given in Eqs. (|A3|) and (|A4[) . The system described by Eq. (|A9|) is equivalent to an Ising 
ferromagnet in an external magnetic field. Thus the point A corresponds to the second order critical point of the 
equivalent magnetic system at 7* = a* = 0. The location of A in terms of the original coordinates is: 

A-^ + ^-^V (All) 

thus A is always located in the region a > and a < 0, as shown in Fig.O The first order line AB in the COE phase, 
determined by the condition 7* = 0, 

AB -- a = -kb a+ ^ + ^ (A12) 

is straight with negative slope, as shown in Fig. O For negative a, r jumps discontinuously across AB between the 
solutions in the COE phase, r = ±y/—a*/b. For positive a, there is a jump between the COE phase and the NOR 
phase as shown below. Thus the coordinates of the point B are 

B-(£ + £..). 

The line BC, obtained from fj( COE ) = fj( NOR ) = at the potential minimum, is first order, since the bracket in 
Eq. (|Aip always vanishes in the COE phase and either a or d must be discontinuous from the NOR phase to the COE 
phase. The expression for BC is easy to derive but is too complicated to show here. For sufficiently large a, the COE 
phase turns into the NG phase, as shown below. The boundary CC separating the NG phase and the NOR phase 
is not affected by the coupling 7 since d vanishes in both phases. Then the condition f2( NG ) = f2( NOR ) = at the 
potential minimum implies 

CC : a = ^. (A14) 

• The phase boundaries BB' and CD' 

To determine the boundaries primarily associated with the diquark condensate, it is most useful to study the reduced 
free energy f2( COE )(cr(<i), d), Eq. (|A5|) . When (3* is positive, the system is equivalent to a usual Ising ferromagnet and 
shows a second order transition at a* = 0. One finds two boundaries, 

BB' : a = 0, (A15) 



CD' : a = Uc + y/c 2 - Aba). (A16) 
6 

We have chosen the larger root of 6a 2 — 2^ca + Ar^a — obtained from a* = in Eq. (IA16|) as the only root 
compatible with the condition that l^ 00 ^ = fi( NG ) < fi( NOR ). Thus the coordinates of the critical end point C are 

c = (w 4 ir) for "' >a (A17 » 

The condition (3* > on the line CD' is equivalent to < 7 < (c/Z)y/]jjb. If (3* < 0, or equivalently 7 > (c/3)y/fi/b, 
a tricritical point D determined by a* — (3* — appears on CD', as shown in the right panel of Fig. [2l In this case, 
the line CD becomes first order and C becomes a triple point. The location of D is 

V 86 f3 ' 26 J ' v ; 
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with d = \/ c 2 + 16&7 2 //3. The point C is determined by the coexistence of the three phases, COE, NG, and NOR, 
so that the free energies f2( COE ) = f2( NG ) measured with respect to the normal phase, vanish at their minima; then 



[ 96 'V b \3Vb VP) 



for 0* < 0. 



(A19) 



b. Negative b 

For b < 0, the presence of the sixth order /-term in Eq. (fl7|) is crucial. Nonetheless, the topological structure of 
the phase diagram in Fig. [5] is exactly the same as for b > 0. 

• The critical point A and the first order phase boundary ABCC 

With b < and the / term, the free energy f^ COE )(cr, d(a)) in the COE phase has at most three local minima, 
which we label &i (i — 1,2,3). The competition among the three phases leads to at most three first order lines. 
However, the detailed stability analysis outlined below show that there is only one first order line and critical point 
in the COE phase. Thus the phase structure for b < is the same as that for b > 0. 

Let us now determine the two competing chiral condensates in the COE phase, labelled by a% and a 2 , at which 
n[ COE) = n { 2 COE) = J2 min . Then fi( COE )(M(<j)) = -a 2 /4fj + 7 V + a*a 2 /2 - ca 3 /3 + ba 4 /4 + f<j 6 /6 = (//6)(<r - 
<7i) 2 (<7 — a 2 ) 2 (a — ti)((T — T2) + f^min where t\_ 2 arc complex except when t\ = T2 = 03. Eliminating and defining 
r = a 1 + (T2 and s = c t \<J2 1 we have 



and 



4/r 3 + (3b - 6/s)r - 2c = 0, (A20) 

= fr A + h -r 2 - fs 2 + bs, 

fr 2 ~f S +^j, (A21) 

a„i„ = - f -(2s-3r 2 )s 2 . 

6 V ' 4 4/3 

In the COE phase, several conditions must be fulfilled: (i) a h2 > a/27, (ii) fix = Q 2 < ^3, (hi) ^ (COE) < n (NOR \ 
and (iv) fi( COE ) < fi( NG ). From Eq. (1^201) . we derive a necessary and sufficient condition for (i)-(iv), < 2y/s < r. 
Under this condition, Eq. (|A20|) . a cubic equation in terms of r, has only one root satisfying r x < r < r c . The upper 
limit r x and lower limit r c , and correspondingly s x and s c , satisfy 

5/r 3 + 66r c - 4c = 0, s c = 4r 2 , (A22) 
4/r 3 + 3br x - 2c = 0, s x = 0. (A23) 

The root r = r c corresponds to the critical point A since the two minima merge into one, <t± — a 2 , while r = r x 
corresponds to the point B where the two local minima are a\ = and <j 2 7^ 0. Using Eqs. (|A20[) , (|A21[) and (|A22|) . 
we thus find 



(26r c 3 -3cr c 2 )^), 

B ^ {-\^ + h + 2 -f^- (A25) 

It can be shown that the point A is always in the fourth quadrant (a > and a < 0) and that the point B is always 
on the a-axis. 

The line AB located in a < corresponds to the interval, r x < r < r c . Thus there is only one first order line, AB, 
and critical point, A, in the COE phase. The line BC, whose explicit form is complicated, is first order, as for b > 0, 
since Eq. (|A1[) is always satisfied in the COE phase and either a or d is discontinuous across BC. The boundary CC 
is not affected by the diquark condensate and, following the similar argument leading to Eq. (|A20[) , we have 

CC : a = - h -v 2 + C -r x , (A26) 
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which is in the region a > 0. 

• Phase boundaries BB' and CD' 

As for b > 0, the lines BB' and CD' obtained from a* = are 

BB' : a = 0, (A27) 
CD' : fa 4 + 4j 2 ba 2 - 8 7 3 ca + 16 7 4 a = 0. (A28) 

The maximum root of a in Eq. ([A28]) should be chosen to satisfy the condition on CD': ft( COE ) = fi( NG ) < f^ NOR ). 
The coordinates of C for 8* > are easily found from the crossing of CD' and CC: 

C= f-|r x 2 + |r X) 27r x ), for 8* > 0. (A29) 

If / < 6 3 (c — 3c")/ (c + c") 3 with c" = yj c 2 — 8&7 2 //3, /3* can change sign on the line CD', and a tricritical point D 
determined by a* — 8* = appears, as shown in the right panel of Fig. [5] In this case, the line CD becomes first 
order and the C becomes a triple point. 

Thus we find that the topological phase structure in the massless three-flavor system with 7 > and A = is 
independent of the sign of b, and a critical point A always appears in the COE phase as shown in Fig. [21 



c. Effect of small X 

Finally, we discuss the effect of the term Ad 2 er 2 , hitherto neglected. With this term, the stationarity condition 
becomes [cf. Eq. (|A1[) ], 

dn 3F a 8 2 2 

__ = _ + _ d - 1G + X a =0. (A30) 
Accordingly, the reduced free-energy fl^ COE \a,d(a)) in Eq. (|A2|) is modified as 

n^(a,d(o-)) = + j*o + fa 2 - ~~<7 3 + + {a 6 , (A31) 
where the new coefficients are 

a* x = a*-2a^, c x = c - 67^, 6 A =6-4A-^. (A32) 

Therefore, insofar as \X/8\ <C 1, as is suggested microscopically (see Appendix [Cj , the A term does not qualitatively 
change the phase diagram in Fig. [5] 



2. Two massless flavors 



Following the same argument as in the three flavor case, we use the stationarity condition for the GL free energy 
n 2F in Eq. dHJ>, 

dn 2F OJ ( a , 8 2 2 



dd V 2 2 

to obtain the reduced free-energies: 



2d - + ^-d 2 + Act 2 = 0, (A33) 



for cr 2 < -a/(2X), and 



tt( COE >(M(a)) = - g + ^a 2 + ^ + {a 6 , (A34) 



«( COE )(a(d),d) = fi + ^rf 2 + ^d 4 - |Jd 6 , (A35) 
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for d < —a/ (3. Here the coeficients are 

a' = a — 2Xa/j3, 
b' = b-4\ 2 /{3, 

nl F = a(-12A 2 a + 3A6a-/a 2 )/(48A 3 ), (A36) 
a' = -{fa 2 -2Ak* + 4A 2 a)/?7(8A 3 ), 
& = (b'\-fa)f3 2 /(4\ 2 ). 

Note that the potential (|A35[) is bounded, despite the negative d 6 term, a consequence of the condition d 2 < —a/ (3. 

a. positive b 

As in an Ising ferromagnet, the phase boundaries, between NOR and NG, and NOR and CSC, for b > and / = 
are characterized by second order lines, a = and a = 0, respectively. For both /3' > and b' > 0, or equivalently, 

i^>^>0, (A37) 

the COE-NG and COE-CSC phase boundaries are characterized by a' — and a' = 0, respectively. These conditions 
can be rewritten as 

a = COE - NG boundary, (A38) 

a = ^— \ COE - CSC boundary, (A39) 
2A 

depicted in the left panel of Fig. [3J 

On the other hand, when A > y/U/3/2, the COE phase ceases to exist. Then the competition between the free 
energies = — a 2 /(46) and fij^jn = — a 2 /(4/3) determines the first order NG-CSC boundary: 

a = ay/f3jb, (A40) 

as shown in the right panel of Fig. [3] 

b. negative b 

The case b < requires introduction of a positive / term. 

• First-order phase boundary E'EFF' 

The phase boundaries EE' and FF' are determined by fi( NG) = fi( NOR ) and fi( COE ) = fi( csc ) at the potential 
minimum, 

EE' : a = 36 2 /(16/). (A41) 



FF ' : -s 



a-3(fc-4A 2 //J) 2 /<16/)l . (A42) 



To determine the boundary EF, we use the free-energies of the NG and CSC phases, 



f2 (NG) = 

mm 

^" C) = -"7(4/3), (A43) 



6bfa + (b 2 - 4/a) 3/2 - b 3 j (24/ 2 ) 

where the global minima at a ^ in the NG phase and d ^ in the CSC phase implicitly ensure b 2 — 4/a > 

i 

{/3 6b fa + {b 2 - 4/a) 3 / 2 - b 3 j (6/ 2 )} V2 . (A44) 



(a < 0). From = fij^n > we obtain 



EF : a 
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The critical end points E and F are then 



B_(£, oV F-fVfv + a*),**). (A45) 



\Wf J ' \Wf V f J' 2/ 

• Second order phase boundary EG' and FH' 

The boundary EG' is not affected by the chiral condensate and is characterized by a — 0, while FH' is given by 
a' = 0: 

EG' : a = 0, (A46) 
FH' :a=j(b- y/b 2 - 4/a) . (A47) 

These boundaries indeed terminate at the points E and F given in Eq. (|A45| . On EG' and FH', the coefficient /?' 
never changes sign, and therefore a tricritical point does not appear, unlike in the three-flavor case (the point D in 
the right panel of Fig. 

While we have left out details of the derivation of the phase structure in Fig. [H the structure can, as mentioned 
earlier, be explicitly checked by direct comparison of the free energies. 



APPENDIX B: SIGN OF 7 FOR THREE FLAVORS 



In this Appendix we summarize how the anomaly-induced interaction —jd 2 a in Eq. (JTTJ) results from the basic 
6-Fermi interaction. The latter has the form [3] 



H 



inst 



— — Arist — — .9D 



detr a (l- 75 )g£ + h.c. 



4et QjtaTLa + h.c. 



(Bl) 



with go < 0. The final (co) term of the GL free energy in Eq. ([7]) originates from the expectation value (in lowest 
order) of Hi ns t; inserting $ from Eq.([4]) into Eq.([7]), we have 



"inst = T G X ( det (^.a9lJ + h ' C -) = T ' G X a ■ ^k^yk-({^ Ra 4aW } Rb <li b )(t Rc qi c ^ + C.C.. 



(B2) 



Thus the GL coefficient is related to the microscopic coupling by c G x /2 = — 2gr> > 0. The coefficient Co must be 
positive in order for m 2 , ~ Co to be positive (see Eq. (jDip ). In addition, for cq > 0, the chiral condensate at low 
temperature is positive. We have G x > as well. 

To proceed, we replace one of the q R qL pairs in Eq. (|B2p by the chiral field which gives 



CO 



T ■ G X ' Z ' e Vk£V 3 >k>{{? Ra qia){q R b < lLb))®k>k + C.C. 

3 2 
6 ' 4 



c 3 2 

f ' G x ' a ■ - ■ ^3k^yk>{{? R aCq Rb ){qi b Cqi a ))<S> k , k + c.c. 



CO 



Tr[(d H 4)$] +c.c. eee 7l Tr[(d fl 4) $ ] + c - c - 



(B3) 
(B4) 

(B5) 



In deriving (1B4|) from (|B3j) we use the fermion anti-commutator twice and the Fierz transform; The overall sign 
changes from — to + because the coefficient of the CC term in the Fierz transform is negative. Then in deriving (1B5|) 
from (|B4[) we use the definition of the diquark condensate, Eq. (J5]). The contraction of the e-tensor three-times gives 
+ (— ) 2 preserving the overall sign. The coefficient 71, defined in Eq. (fl4| is thus co(G x /Gd) 2 > 0. 

With the CFL ansatz lfT6|) . the anomaly-induced term reduces to 

^inst — > -37id 2 CT = ~~/d 2 a. 
The minus sign originates from = —d R . Thus we finally find that 7 [Eq. (|17p ] is given by 

rA 1 >0. 



(B6) 



7 = 3 71 =3 C (^ 



(B7) 



This result is consistent with a similar result given in [17j | , where it is found that the effect of the quark mass combined 
with the axial anomaly reduces the energy of the CFL phase. 
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APPENDIX C: A, f3 AND b FROM THE NJL MODEL 



The A in Eq. (fT7|) may be evaluated in the three-flavor and two-flavor NJL models as follows. We start with the 
free energy Eq. (5.15) in Ref. [8] for three flavors, and Eq. (4.50) for two flavors, in the chiral limit: 



^3f(ct, d) 



3a 2 
2G, 



co 



2£r v 2G,i 



3 3d 2 f d 6 p ^ r 

4w ± 



^ I d 3 p ^ r 
2G d J (2^)3^4 



8^5 



2E* 



f + 16Tln (l + e^ /T ) + 2T\n (l + e~^ /T ^ ,(C1) 



8T In 1 



(l + e - w± / T )+4Tln( 



4Tln (l + e- E± / T 



where 



E ± = ^TW 2 ±p with M N-coG x a 2 /2 (three flavors) 
v r n { a v two fl avors ) 



± 



d 2 

U 8(l) 



UJ 



±\2 



d 2 . 



(C2) 

(C3) 

(C4) 
(C5) 



Here = d is the octet gap, and d\ = 2dg the singlet gap, for nine quasi-quarks (three colors times three 
flavors). The Co-term originates from the instanton- induced interaction; we ignore in Eq. (|C1[) the contribution from 
the quark— quark and quark— anti-quark six-point instanton-induced interaction (the 71-term in Eq. (|14|) ). Since the 
fcrmion integral does not diverge in the infrared, a naive power-series expansion of f2 in d and a should be valid. 

The gap is determined by minimizing il with respect to d 2 : 



3F 



dd 2 

0Q.2F 

dd 2 



3 
2G~d 
1 

2G^ 



E 

± 

E 



d 3 p 
d 3 p 



ui, 



4 , ^8 

-rr tanh -2- 
± 2T 



2 , Ujf 

— - tanh — — 
2T 



2 . 
— r tanh 



± 



2T 



(C6) 
(C7) 



where we introduce an ultraviolet cutoff at p = A to regulate the integrals. In the NJL model, A ~ lGcV, while in weak- 
coupling QCD at high density, the momentum dependence of the gap function leads to A ~ 2567r 4 (2/3g 2 ) 5 / 2 /i ^s> [1 
[l5T |. In the following, we always assume A>(i and A » T. The — terms are the contribution from the particles 
and the + terms from anti-particles. The critical temperature T c of the super-to-normal transition for three flavors is 
determined by setting the right side of Eq. (|C6[) . and of Eq. (|C7I) for two flavors, to zero. 
The term Xa 2 d 2 in the GL expansion of Cl(a, d) is obtained as 



A3F — 



'3F 



da 2 dd 2 



a=d=a 



—T 

2ir 2 V 







dpp^r- 
Op 



1 



p±H 



■ tanh 



(p±aO 

2T 



(C8) 



Evaluating the above integral by partial integration, making the approximation, tanhx ~ 1 (x > 1), tanh x ~ — 1 [x < 
— 1) and tanhrr ~ x (—1 < x < 1), we arrive at 



\3F 



2tt 2 



hi 



A 
2T 



> 0, A 



2F 



A 3 f/3. 



We can similarly calculate the coefficient (3 in the term (3(d ) / '4. For three flavors 



/3 3 f 



a 2 r2 



3F 



ad 2 ad 2 



G 



cr=d=0 



which leads to 



@3F ^ -o 



f — 

V2T 



dP 



2T 



9 



p± fj, dp 



p±H 



tanh 



(p±/*) 
2T 



>0, /3 2F = /3 3 f/6. 



(C9) 



(CIO) 



(Cll) 



For typical values in the NJL model at intermediate density: A = 1 GeV, p = 500 MeV and T = 50 MeV, we 
have A3F//33F — 0.7/18 = 0.04. Also, at asymptotically high density in QCD with A cx p, we have \sf / (3tsf ~ 
ln(/x/T) 2 / (p/T) 2 which is parametrically small at large /i. 
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The coefficient b in the term b(a 2 ) 2 /4 can be calculated from the free energy Eq. (|C1[) . or Eq. (|C2p , with d = 0: 



6„. = 2 - - - 



da 2 da 2 



3N f 



,A 

_i_ Jo 



dp 



1 , (p±m) 

- tanh — — — 
p 2T 



(C12) 



a=0 ± 

where we have assumed M = a for both two and three flavors, for simplicity. Evaluating the integral, we find 



3iV 



4tt 2 



In 



A/2T 



1 + /i/2T 



i-JL) 



In 



A/2T 



l-/i/2T 



(C13) 



Unlike for the coefficients A and /3, the sign of b changes from positive to negative as \x increases from zero. 

Finally, we show the integral form of the coefficient b 2F which appears in the term, b 2F d 2 a i /4 in the GL expansion 
of the two flavor NJL model: 



d 3 il 



b2F 2 da 2 da 2 dd 2 



a=d=0 



ld_ 

p dp 



1 tanh^ 



p±H 



2T 



(C14) 



where the sign of b 2F changes, as one can show, from negative to positive as y, increases from zero. 



APPENDIX D: MASS SPECTRA OF THE 77' MESON 




FIG. 7: Mass spectra of 77' meson in the intermediate density region when (a) a-d mixing does not exist (7 = A = 0); (b) 
a-d mixing is considered (7 / 0,A / 0); (c) a finite quark mass is considered (m q 7^ 0) in addition. Here m^, = 36cocr 3 //y, 

m \ 2 ( m <? = 0) = 12 (xi q= V x i - and 2 = 12 ( :E3 T (^i^s - x^)/y/x\ - X2J m q . 



The 77' mass spectra can be found, similarly to that of the pions, by diagonalizing the mass matrix Eq. ((5T|) . Figure[7| 
shows the mass spectra for 77': (a) with neither the a-d coupling nor a quark mass; (b) with the a-d coupling but 
without the quark mass; and (c) with both the a-d coupling and non-zero m q . In case (a), fj is massless, since there is 
no term that breaks the U(1)a symmetry in Eq. (jT3j ; on the other hand, is massive because of the axial anomaly, 
with 

a 3 

m 2 , = 36co^2". (Dl) 

J-q' 

This result differs notably from that for pions, where both 7r and tt are massless in the absence of the a-d coupling 
and the quark mass. 

With the a-d coupling alone (b), the coupling mixes 7/ with fj' and makes the mass eigenstates r)' x 2 much heavier 
than in the absence of the coupling. We have 

™< 2 K = 0) = 12 [x x T y/xl-x 2 ) , (D2) 
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where 

2 7id 2 cr + \ 3 d 2 a 2 A^d 2 a + X 3 d 2 a 2 
X! = m„, + - 2 + 3 , 

J rj' J fj' 

Ac 0ll d 2 <T A + c \ 3 d 2 <7 5 + 6 7 iA 3 d 4 a 3 
^2 = 6 -2—2 . (D3) 

The 71-term not only mixes the mass spectra of rf and ff , but also serves itself as a mass term for 77' and ff. On the 
other hand, the A3 term, which preserves U(1)a symmetry, acts only to mix rf and ff. 
In case (c), the non-zero quark mass increases both the rj'i and rf 2 masses by 



m 



ml , (m, = 0) + 12 Us T t^J* m (D4) 



with 



■I'3 



A cr 4rid 2 
— - — 1 — - — 

f'2 ' f 2 1 
J rj' J fj' 

6c rid 2 (7 3 + 4rid 2 ( 7 id 2 cr + A 3 d 2 ff 2 ) + A cr (4 7 id 2 a + A 3 d 2 <J 



2\ 



* 4 = 2— ^ ^—L "—^ L (D5) 

For small m 9 , ff x becomes heavy owing to the a-d coupling. By contrast, the pion 7ri remains light, even in the 
presence of the coupling. 
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